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Abstract 



Systems of Newton equations of the form q = — 2 A~^(g)V/c with an integral 
of motion quadratic in velocities are studied. These equations generalize 
the potential case (when A=I, the identity matrix) and they admit a curi- 
ous quasi-Lagrangian formulation which differs from the standard Lagrange 
equations by the plus sign between terms. A theory of such quasi-Lagrangian 
Newton (qLN) systems having two functionally independent integrals of mo- 
tion is developed with focus on two-dimensional systems. Such systems 
admit a bi-Hamiltonian formulation and are proved to be completely in- 
tegrable by embedding into five-dimensional integrable systems. They are 
characterized by a linear, second-order PDE which we call the fundamental 
equation. Fundamental equations are classified through linear pencils of ma- 
trices associated with qLN systems. The theory is illustrated by two classes 
of systems: separable potential systems and driven systems. New separation 
variables for driven systems are found. These variables are based on sets of 
non-confocal conies. An effective criterion for existence of a qLN formula- 
tion of a given system is formulated and applied to dynamical systems of 
the Henon-Heiles type. 



I. Introduction 

In this paper we introduce and study such systems of Newton equations 
q = M{q) that can be generated as equations of the form 

^_ d dE dE _ 
dx dq dq 

by an energy-like function quadratic in q 

n 

Eiq, q)= J2 ^iM)m + Ki) = (I'M + Kq). (1-2) 

where A{q) is an n x n symmetric matrix with real entries Aij{q). Here and 
in what follows we use the standard mechanical notation q = (gi, . . . , 
q = ((71, ... , (?n)*, for position and velocity vectors (the superscript t denotes 
the transpose of a matrix), where qk = -§^qkT k = \, . . . ,n with x G R being 
the independent (time) variable. By Newton equations we mean second 
order ordinary differential equations (ODE's) of the form: acceleration q is 
equal to the velocity independent force M{q). The force M may be potential 
or not. 



The equations in (|l.l| ) are called here quasi-Lagrangian (qL) equations since 
they differ from the Lagrange equations for E{q,q) by sign between terms 
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only. These equations are shortly denoted 
where 



= {6tE,...,6tEY 



d dE 
dx dqi 



+ 



dE 
dqk 



The qL equations are not invariant with respect to arbitrary point transfor- 
mation but it can be easily shown that they remain invariant with respect 



are 



to the affine change of variables q = SQ + h where Q = {Qi, . . . , Qn] 
the new variables and S £ GL{n), h G R". 

In the present article we shall mainly discuss quasi-Lagrangian sets of New- 
ton equations (qLN) generated by a function E of the form ( |1.2| ) in the 
two-dimensional space of variables q = {qi,q2) = {r,w). This class of equa- 
tions (which seems to be completely new) is a very interesting class because 
of its rich differential-algebraic structure and also because it contains (as 
special cases) the well understood class of point-separable potential Newton 
equations q = —dV{q)/dq and the class of non-potential Newton equations 
of the triangular form r = Mi(r, lu), iv = M2{w) which we shall call driven 
systems. The qLN systems are not necessarily Lagrangian and thus they do 
not have any straightforward Hamiltonian formulation. 

In this paper we develop a theory of completely integrable sets of qLN equa- 
tions characterized by the existence of two functionally independent integrals 
of motion quadratic in velocities: E as above and F = q*B{q)q + l{q). The 
existence of a second integral of motion has far-reaching consequences; it 
eventually leads to wide classes of completely integrable qLN systems. 



Example 1.1 The function E 
inserted into (|l.l| ) gives rise to 



rrw 



wr 



awr"^ _l_ Ifir^ + ^ when 



d dE _|_ dE_ 
dx dr ' dr 

d dE I dE 
dx diii dw 



-2w{r 



-2w r 
r 



ar + ^) -|- r{w — Aaw + d) 
r{f -ar + ^) 



r — Ml (r, w) 
w — M2{w) 



which is equivalent to a set of two Newton equations 



r = 
w 



ar — 
- 4ati; 



: = Mi{r,w) 
d = M2{w) 



since the matrix 



-2w r 
r 



(1.3) 



(1.4) 



is nonsingular. We see that the operation Q = 5^E generates linear combi- 



nations of the Newton equations (|1.4D . 
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Equations ( |1.4| ) were discovered accidentally as a Newton parameterization 
of the second stationary flow of the Harry Dym hierarchy Q: 




(here d = d/dx) where we substituted u = r~^. The substitution w = 
—r^r + ar^ gives the system ( |1.4D . The particular feature of (|1.4| ) is that 
it is a driven system: the equation for w can be solved independently and 
then the solution w{x) drives the equation for r. 



II. General properties of quasi-Lagrangian Newton 
systems 

Let us consider an n-dimensional qL system = 5^E with (quadratic in 
velocities) energy-like function 

n 

Eiq, q)=J2 ^iji<l)mj + Hq) (2.1) 

with a symmetric (which can be assumed without loss of generality) matrix 
A(q) = A^{q). We shall formulate the necessary and sufficient condition for 
the matrix A{q) to make the equations = 6~^E equivalent to the set of 
equations 

= ^' - M{q) (2.2) 
with a velocity independent force M{q) = {Mi{q), . . . ,Mn{q)Y- 



Theorem 2.1 For the function E given by (2J_) with a nonsingular matrix 
A{q) the following conditions are equivalent: 

1. The equations ^ = 5^E are equivalent to the set of Newton equations 
q = M{q) with velocity independent forces M = —^A~^{q)Vk{q). 

2. The function E is an integral of motion for the qL system {) = 5^E 

3. The matrix elements Aij{q) satisfy the following set of "cyclic" diffe- 
rential equations 

= diAjkiq) + djAkiiq) + dkAij{q) for alH, j, A; = 1, . . . , n. (2.3) 
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Throughout the whole article the symbol V denotes the gradient operator 
and di = d/dqi. Later on we will also use the notation dij = d'^ jdqidqj. 

Statement 2 of the above theorem explains the name "energy-like" for the 
function E. 



Proof: Let us calculate the i-th equation in = 8^E: 



d BE dE d 



2 AiM)<ii + d^k + Y, {diAjkiq) + djAki{q) + dkAij{q)) qjqk- (2.4) 



The last equality in ( |2.4D is due to the symmetry of A{q). Thus, clearly, 
2Aq + Vk = if and only if the equations (^^) are satisfied and the equiv- 
alence of 1 and 3 is established. 

Let us now calculate the total derivative of E with respect to x. 



-E' = E I 2 E Aijqj + dikjqi + Y^ dkAijqiqjqk 



E 



2 Aijiij + dik^ qi + ^Y i^iAjk + djA^i + dkAij) m^qk- (2.5) 



The second term on the right hand side of the above equation has been 
rewritten by renaming indices. It contains precisely the cyclic conditions 



(2.3). So, if one (and thus both) of the statements 1 and 3 are satisfied, 
then both terms in (|2.5| ) vanish. On the other hand, if = then terms at 
different powers of qi in (p.5D must be equal to zero, which implies both 1 



and 3. QED 



Remark 2.2 For n=2 the general solution of equations ( |2.3| ) can easily be 
found. It is 

Aii{w) = anP' + bw + a 
2Ai2(r, w) = —2arw — br — cw + (3 (2-6) 
A22ii") = ar'^ + cr + J 

with some real constants a, b, c, a, /3, 7. The corresponding qLN equations 
read explicitly as 



J^dE I dE 

dx dr dr 

J^dE , dE 

dx dw ' dw 



An A12 
A12 A22 



r — Ml (r, w) 
iv — M2{r, w) 
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where 



, ^ / ^ 1 i A dk , dk 



The remaining part of this work is mostly devoted to the case when qLN 
equations q = —^A^^Vk generated by E admit a second (quadratic in 
velocities) integral of motion F{q, q) = J27,j=i Bij{q)qiqj + l{q) = q^B{q)q + 
l{q) which is linearly, and therefore functionally, independent of E. 

Theorem 2.3 (qLN systems with two integrals) Let the qLN system 
of Newton equations 

= 5+E = 2A{q + ]^A'^Vk), (2.7) 

generated by the function E{q,q) = q^A(q)q + k{q), admit a second, func- 
tionally independent quadratic integral of motion F(q,q) = q^B{q)q + l{q). 
Then 

1. The matrix B{q) has the same structure as the matrix A{q) in the sense 
that the coefficients Bij{q) of B{q) satisfy the set of cyclic differential 
equations (^.Sj). 

2. Ifdet{B) / then 

A-^Vk = B-^Vl (2.8) 

and so the qLN system = 5^F = 2B{q + ^B^^Vl) generates the 
same Newton equations as E. 

3. Any differentiable function f{E, F) generates the same system of New- 
ton equations (by = 5~^f{E,F)) as E does. In particular, any linear 
combination \E + fiF generates the same system of Newton equations. 

The statement 2 shows one of peculiar features of qLN systems: all quadratic 
(in velocities) integrals of motion of a qLN system generate the same system 
(see also sec. 151 in |§). 

Proof: The requirement B = yields (cf. (|2.5|)) 

^ = X! p X! + dil\ (ji + Y^ dkBijqiqjqk = 
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= J2'ii (2^ (-^^~'^^) + . + E dkBiMAk (2.9) 

where the index i at the vector expression containing matrices B and 
denotes its i-th component. The equahty is satisfied identicahy with respect 
to q and so both sums must be separately equal to zero. It follows, that the 
Bij satisfy the cyclic conditions diBjk+cycl. = and that 2B (—■^A^^Vk] + 



V/ = 0. The latter yields precisely the equation (|2.8|) since we assumed 
det(i?) 7^ 0. So the statements 1 and 2 are proved. 

The operator 5^ acts as differentiation on the algebra of constants of motion, 
so that 

= 6-^f{E,F) = %6^E+%5+F = 



dE dF 

(q-M) 



\dE dF 

(where M = —^A^^Vk = —^B^^Vl) which proves the statement 3 of the 
theorem. QED 



It is important to stress that the equation ( |2.8| ) is the necessary and sufficient 
condition for the equivalence of the qLN system ( p.Tl) and the qLN system 
generated by F = q*'B{q)q + l{q). This condition will be used later. 



III. qLN equations in two dimensions 

We shall from now on restrict our considerations to the case n = 2. We will 
use the notation q = (^1,^2)* = {r, wY . The case of arbitrary n is studied in 
a separate paper [0]. 



For n = 2 Theorem 2.3 contains two special cases which explain the connec- 
tion of our theory with classical results about separable potential Newton 
equations and with the class of driven systems where one of the Newton 
equations depends only on a single variable r or w and can be solved on its 
own. 

Corollary 3.1 Assume that the Newton equations 

f = Mi(r,w) , w = M2{r,w) (3-1) 

generated by the integral E = q*A{q)q + k{q) (with the matrix A given by 
^)) asO = 6+E have a potential force: Mi = -dV/dr, M2 = -dV/dw. 
Then the potential V{r, w) satisfies the Bertrand-Darboux equation 

= {Vyjw — Vrr){—'2.arw — br — cw + /3) + 2Vrw{a'w'^ — ar'^ + bw+ , . 
-cr + a - 7) + 3Vr{2aw + b)- 3I4,(2or + c) ^ ' 
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(where the indices at V denote partial derivatives with respect to r and w) 
with the coefficients a, b, c, a, /?, 7 being exactly the coefficients of the poly- 
nomials in entries of the matrix A as given by pT^. This means that the 



Newton system ( p. ]\) can be solved by separating variables in the related 



Hamilton- Jacobi equation (see 1^). 



Proof: If M is potential, then according to Theorem 2.3 M = —\A ^Vk 



—'VV and so VA; = 2A'W. The potential V exists provided that d'^k/drdw = 
d'^k/dwdr. This yields exactly the Bertrand-Darboux equation (|3.2| ) for V . 
QED 



Remark 3.2 The quantity k{r, w) / det{A) satisfies the same Bertrand-Darboux 
equation as the potential V. This result can be verified directly but it also 
follows from Theorem |4.1| in the next section. 



Remark 3.3 Let us emphasize that the Hamiltonian system 

dV dV 
r = s, w = z , s = — — , z = - — 

or ow 

generated by a separable natural Hamiltonian H = ^(s^ + z"^) + V{r, w) can 
be reconstructed as the qLN system = 6~^E = 2A{q + ^A~^Vk) from its 
second integral of motion E. This is easy to see, since the above Hamilton 
equations are equivalent r = iu = 

The second class of equations satisfying the assumptions of Theorem |2.3| is 
the class of qLN systems of the form 

r = Mi{r,w) , w = Nhiw) (3.3) 

which naturally generalizes the system in Example Such systems are 
called driven since the equation for w can be solved independently and then 
w{x) can be substituted into the equation for r. Observe that the second 
equation (and thus the whole system) admits an extra integral of motion of 
the form F{w,w) = jl — J M2{w)dw. The qLN system = 5^ E attains 
the form ( p.3p if and only if the second component M2 of the force — ^A~^V/c 
does not depend on r 

|.(^-iVA;)2 = 0. (3.4) 



Example 3.4 The qLN equations generated by the function 

E = rrvj — wf"^ + k{r, w) 
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are driven (i.e. have the form ( |3.3| )) provided that k{r,w) satisfies the fol- 
lowing second order PDE 

a /I 2w 



dr \r 



which is a specialization of (3.4). The general solution of the above equation 
is 



k{r,w) = f 



w 



+ r g{w) 



with arbitrary twice differentiable functions / and g. The corresponding 
qLN system attains the form 

T ( d 

f = -rg'{w) + -^f — , w) = -2— {wg{w)) 



w 



dw 



and can be solved by quadratures (see Section 7). The second integral of 
motion of our system, F = /2 — J M2{w)dw = + 2wg{w), yields the 
matrix B 


1/2 



B 



which is singular so F does not generate our system. However any linear 
combination \E + ^iF of E and F (with both A and /i 7^ 0) is another 
integral of motion with a non-singular matrix B' = \A + jiB generates the 
same driven system as E. 



Existence of two functionally independent constants of motion does not au- 
tomatically imply Liouville integrability since we also need a Hamiltonian 
formulation for our equations of motion. Our systems usually do not have a 
Lagrangian formulation and so they do not have the standard Hamiltonian 
formulation. On the other hand the special system discussed in Example 
1.1 , being a stationary flow of the Harry Dym hierarchy, is expected to be 
integrable. The question thus arises if/when our qLN systems possess a 
non-standard Hamiltonian formulation. In Section 6 we shall demonstrate 
the existence of new Poisson structures for qLN systems and their close rela- 
tionship with Poisson pencils for separable potentials. We shall also explain 
there when and in what sense our qLN systems are integrable. 



IV. Fundamental equation 

We shall now characterize those two-dimensional qLN systems which admit 
two (quadratic in velocities) functionally independent integrals of motion E 
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and F, with the force M 



i-B ^V/. We remind the reader 



that for n = 2 we use the notation q = (gi, (?2)* = {r, u;)*. 
Let us consider two symmetric 2x2 matrices A^rjw) and B{r,w) both 



satisfying the cychc conditions (2.3). According to Remark 2.2 they must 
have the fohowing structure 



" Au 


Al2 


, B = 


' Bn 


B12 


_ Ai2 


A22 


B12 


B22 



A 



with the polynomial entries given by (cf. ( |2.6D ) 

Aii{w) = ai'uP' + hiw + ai 

2Ai2{r, w) = —2airw — bir — ciw + /3i 

A22{r) = air^ + cir + 71 

Bii{w) = a2W^ + h2W + a2 
2Bi2{r, w) = -2a2rw - b2r - C2W + /32 
B22{r) = a2r^ + C2r + 72 



and 



with some arbitrary real constants ai, . 



(4.1) 



(4.2) 



(4.3) 



,72- 



Theorem 4.1 (fundamental equation) Let 



r 
w 



-A-^Vk 
2 



1b- 

2 



(4.4) 



with nonsingular 2x2 matrices A, B given by ^.\), ^■^ ) and l{4-^ , be a set 
of qLN equations. Then the functions Ki = k/ det(74) and K2 = 1/ det{B) 
both satisfy the same linear, second order, partial differential equation 



= 2(^112^22 - ^22^12) i^rr 

— 2(^iii322 — ^22-Bii) Krw 

+ 2{AiiBi2 — A12B11) Kujy^ 

+ 3(^129^^22 - Bi2drA22 + ^229^^11 

- 3(^119,^22 - Bll9,^22 + ^129^^11 
+ 3(9^.^22 d-a^Bii - drB22 9^^ii) K 



B22dwAu)Kr 

-Bi29,i,^ii) Ku 



(4.5) 
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which explicitly reads 



= 2i^'rr[72/5i - 71/^2 + (^271 - 72^1 + PiC2 - cil32)r + (71 C2 - 72Ci)u; 
+ ip2Ci - C261 + a2/3i - ai/32)r^ + 2(7102 - ^2ai)wr 
+ (0162 - a2bi)r^ + (a2Ci - C2ai)z«r^)] 
+4Krw[a27i - "172 + (a2Ci - aiC2)r + (6271 - 72^i)t« 

+(0201 - aia2)r'^ + (71 «2 - l2ai)'w'^ + (hci - C261)™ 
+ (ai^2 - a2bi)wr'^ + (a2Ci - C2ai)rt(;2] 
+2iC^^„[ai/32 - 02/?! + (0261 - 0162)?' + (a2Ci - aiC2+ 
6i/32 — b2Pi)w + (ai/32 — 02/^1 + fe2Ci — C2bi)uP' 
+2(0201 - aia2)tfr + (0201 - C2ai)w'^ + (0162 - a2bi)rnP'] 
+3i^^ [26271 - 27261 + /3iC2 - ci/32 + (362C1 - 3C261 + 2a2/?i - 2ai/32)r 

+4(7102 - 72ai)tD + 4(0162 - 0261)^2 + 4(o2Ci - C20i)r7i;] 
+3K^[2a2Ci - 2aiC2 + 6i/?2 - 62A + (2oi/32 - 2o2/3i + 362C1 - 3c26i)w 
+ 4(a20i - ai02)r + 4(o2Ci - C2ai)w'^ + 4(0162 - 026i)r?i;] 
+6iC[62Ci - C261 + 2(0162 - 026i)r + 2(o2Ci - C20i)w] 

(4.6) 

with K denoting either Ki or K2 and Kj. = dK/dr, K^r = d'^K/dr'^ and so 
on. 



Conversely, any solution K2{q) of the equation (4-5) generates two dif- 
ferent systems of qLN equations q = —^A~^\7ki = —^B~^\7h and q = 
— ^A^^Vk2 = —^B^^'Vl2, where the functions ki,k2,h,l2 o,re determined 
by the equations 



h = K2 dei{B) Vki = AB-^V{K2 det{B)) 
k2 = K2 det(A) V/2 = BA-^V{K2 det{A)) 



-lw.^^_.„...^^ (4-7) 



We will call the equation ([4. 51) the fundamental equation associated with the 
matrices A and B. 

The fundamental equation plays a crucial role in our theory of qLN systems. 
Observe that it is invariant with respect to the transformation A 1— > XA+fiB, 
B ^ \' A + {X, X' , fi, fi' G R) since the coefficients at every monomial 
in this equation are skew symmetric in A and B. This is consistent with 



statement 3 of Theorem 2^ , which asserts that if any pair E, F of functions 
generates a qLN system then the linear combinations XE + fiF and X'E + 
fi'F also generate the same system. This explains that the assumption of 
nonsingularity for both A and B is nonessential since if det(^) 7^ a singular 
matrix B can always be substituted by an invertible matrix B' = XA + ^B. 
We shall investigate further properties of the fundamental equation in the 
next theorem. 



Notice that in the second part of Theorem 4.1 one has to reconstruct I2 



and ki by integrating the expressions for V/2 and for V/ci. This can always 
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be done, as the above theorem imphcitly states. Also, notice that in the 
fundamental equation (4^) all terms of degree 4 and higher cancel so that 
the polynomial degree of coefficients in this equation is less than or equal to 
3. 



Proof: (of Theorem |4.lD Our qLN system ( [4.4D is generated by either of 
the two functions E{q, q) = q*'Aq + k and F{q, q) = q^Bq + / and so the 



condition (2^), i.e. A~^\/k = B~^Vl, must be satisfied. This implies that 
V/ = BA"^\/k. This equation for the function I has solutions if and only if 
its compatibility condition Ij-w = Iwr is satisfied. This yields a PDE for the 
function k which, after the substitution k = Ki det(^) and with use of the 
cyclic conditions (|2.3| ), yields that Ki satisfies equation (|4.5| ). By inserting 
into this equation the explicit form of the polynomials An, . . . , B22 we obtain 
(4.6). On the other hand, the condition (|2.8D implies also VA; = AB^^Vl^ 



and its compatibility condition kj.yj = k^j. gives a PDE which in terms of 
K2 = 1/ det(i?) must attain the form ( [4.5| ) with interchanged entries of A 
and B (since the equation VA; = AB~^Vl becomes V/ = BA~^Vk when 
one exchanges A,k and B,l). Due to the skew-symmetry of coefficients of 
the equation for Ki with respect to the entries of matrices A, B (clearly 
seen from the form of ( |4.5D ) the obtained equation for K2 differs from the 
equation for Ki by a minus sign on the right-hand side only. This proves 
that Ki and K2 both satisfy ( ^ ) (notice, however, that this does not imply 
Ki = K2). 

The existence of ki (i.e. the possibility of integrating the equations ( [4.7] ) in 
order to obtain ki) follows from the fact, that the condition d'^ki/drdw = 
d'^ki/dwdl together with V/ci = AB^^V {K2 dei{B)) yields precisely the 
fundamental equation for K2 which is satisfied due to assumptions. One can 
similarly prove the existence of I2 ■ The second statement of the theorem can 
now be proved by checking that both pairs ki,li and k2,l2 given by ( |4.7| ) 
satisfy the condition (|2.8| ) and thus give rise to two systems of qLN equations. 
QED 



Remark 4.2 For B{q) = i/ (a 2 x 2 identity matrix) the equation (4.5) 



becomes the Bertrand-Darboux equation (3^) characterizing all separable 



potentials since in this case q = —\B ^Vl = —VI is a potential equation. 



The next theorem shows that there exists a recursive relation between two 
different qLN systems constructed from a given solution K2 (q) of the funda- 
mental equation ( |4.5| ). This makes it possible to construct a doubly infinite 
sequence of qLN systems corresponding to a given fundamental equation. 

Theorem 4.3 (recursion theorem) Let ki,li and k2,l2 be two pairs of 
functions determined by a given solution K2 of the fundamental equation 
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(4-5) as in ( [^. T[ ). Then these functions are related by the following linear 



algebraic equations 

k2 = h det{AB-^) , h = hTi^AB-^) - h (4.8) 
(where Tr denotes trace of matrix). Moreover, in the infinite sequence 

ko ki k2 

/ / / 

••• i Ki i K2 i K-i ■■■ (4.9) 

/ / / 

Iq li I2 

of triples {Km, km,lm), € Z defined recursively by 

km = lm-1 det(AS-i) , Im = lm-iTT{AB-'^) - k^-l (4.10) 

and by 

Km = km/det{A) = lm-i/det{B) 

the functions km and Im satisfy A^^Vkm = B^^Vlm and thus they both de- 
termine the same (for a givenm) qLN systemq = —^A~^\/km = —\B^^\/lm- 
All functions Km satisfy the fundamental equation ^.3^ ) and are related 
through the following two-step recursion 

Km+i = KmTr{AB-^) - Km^i det{AB-^). (4.11) 
The above recursion is reversible. The solution Km placed between lm~i and 



km determines both Im-i and km- The recursion (|4.1l| ) is soluble. Namely 
if we denote the eigenvalues of the matrix AB~^ by Ai and A2 then it can 
be proved that for the case Ai / A2 the solution of (4.11) is 



Km = T ^{Ki - X2Ko)XT + T ^(i^oAi - i^i)A^ 

Ai — A2 Ai — A2 

while in the case Ai = A2 the solution of ( [4.11| ) becomes 



Km = KoXT + - UA™ 



Ki 
Ai 



1 • 



In both cases Kq and Ki are two subsequent solutions of the fundamental 
equation in the sequence p.9|) which are related by 



V {Ki det{B)) = BA-^V {Kq det{A)) . 
In order to prove the recursion theorem we need the following lemma. 



13 



Lemma 4.4 Let X = AB ^ with matrices A, B as above. Then 

X-^V{det{X)) = V(Tr(X)) 

This lemma follows from the cyclic properties (^j^) of matrices A and B by 
a lengthy but straightforward calculation. 

Proof: (of the recursion theorem) Consider a solution K2 of the fundamen- 
tal equation and the functions ki,li;k2, h defined by (|4.7|). Then obviously 
^2/ det(^) = h/ det(i?) which immediately implies k2 = h det(AB~^). Let 
X = AB-\ Then 

Vl2-VihTr{X) -ki) = 

= X-^V{K2 det(A)) - V{K2 det{B)Tt{X)) + XV{K2 det(5)) = 
= X-^V{K2 det{A)) - {Tt{X)I - X)V{K2 det{B)) - K2 det(S)V(T¥(X)) : 
= X-^V{K2 det(A)) - det(X)V(i^2 det(S)) - K2 det(5)V(Tr(X)) = 
= K2dei{B) {X'^V {dQt{X)) - V(Tr(X))) = 

where we used that -Tr(X)X + det(X)/ = as follows from the Cayley- 



Hamilton theorem. The last equality is due to Lemma [4.4| above. Thus 
I2 = IiTt{X) — ki up to a non-essential additive constant. This proves the 
first assertion of the theorem. 



If we now define the sequence {{km, Im)} via the recursive procedure ( [4.10 ) 



then a simple induction argument shows that each pair (^rrii^m) satisfies 



the condition (2.8) and thus both km and Im determine the same qLN sys- 
tem. Moreover, each Km = km/ det(A) = Im-i/ det(i?) is a solution of the 
fundamental equation as theorem ( |4.1| ) states. Finally, to obtain ( |4.11| ) it 
is enough to insert the formula Km = km/ det{A) = Im-i/ det(i?) into the 
second equation in (4.10). QED 



Example 4.5 (cf. Remark 
sion (14.81) takes the form 



For B = (the potential case) the recur- 



k2 = Wi det{A) , V2 = 2Tt{A)Vi - ki 



with Vi = h- This is the separable case when ( f4.5[ ) reduces to the Bertrand- 
Darboux equation. In the generic case, i.e. when a 7^ in ( |2.(]| ) the matrix 
A{q) can be reduced (with the use of affine transformations q = SQ + h with 
5 G GL(2,R), h e R^, see also Section 5) to the form 



Q1Q2 



9192 

-91 + Ai 



If we now start with the harmonic oscillator potential Vi = ^{qf + 92) then 
the condition Wi = ^A'^Vki gives ki = A291 + ^loi ™d the recursion 
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formulas specify to 



k2 = 2{ql + ql){\i\2 - Aagf - Xiql) 
V2 = Wi + A292 - (9i + <i? 

thus reproducing the potential of the Garnier system Q. It can be shown 
that the above formulas prolongate to the n = 2 case of the recursion for 
the Jacobi family of elliptic separable potentials [^. 



In order to explain the character of the recursion (^|^) more completely let us 
consider instead of the pair (A, B) of cyclic matrices another pair {A+^iB, B) 
with G R. As it can be shown (see below) this pair determines the same 
fundamental equation as the pair [A^B) does. By choosing a solution K2 of 
the fundamental equation and the pair {A + ^iB, B) we arrive at a different 
qLN system q = M^{q) = -\{A + fiBy^\/{K2 det{A + fj.B)). It turns out 
that the force Mfj_{q) is a linear combination of two neighboring forces in the 
sequence (4.9) generated by K2. 



Lemma 4.6 Let A and B be two 2x2 matrices satisfying the cyclic con- 
ditions l \2.3l ) and let K be a solution of the fundamental equation associated 
with A and B. Let also G R. Then 

{A + fiB)-^ V {K det{A + nB)) = A'^V {K det{A)) + /xB"^ V {K det{B)) 



This lemma is a consequence of Lemma iA. It says that a solution of a 
given fundamental equation determines the force M (and so the system of 
qLN equations) up to linear combinations of two consecutive systems in the 
recursion ([4.91) . 

As we have mentioned the matrices A and B uniquely determine the fun- 
damental equation. The choice of A, B which generate a given fundamental 
equation is however not unique since the pair A' = a A + (3B, B' = ^A + 6B 
determines the same equation. One can also ask, to what extent a given 
fundamental equation determines the pair {A,B). The precise relationship 
between pairs {A, B) and the fundamental equation is explained in the fol- 
lowing theorem. 



Theorem 4.7 Let {A, B) 



be a pair of linearly independent matrices A, B 
Then there is a 1-1 relationship be- 



satisfying the cyclic conditions (2.S). 
tween the linear span {XA + ^B : A, /i G R} of A and B and the fundamental 
equation ((.1) i.e. 



1. Any two linearly independent matrices A' = a A -\- I3B, B' = 7 A -|- 5B 
determine the same fundamental equation as {A, B) does. 
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2. If the pair (A' , B') determines the same fundamental equation as {A, B) 
does, then the matrices A' and B' belong to the linear span {\A + ^iB} 
of A and B. 

Proof: An easy calculation shows that the fundamental equation associated 
with the matrices A' = aA + PB and B' = 'jA + 5B differs from the funda- 
mental equation associated with the matrices A and B by the multiplicative 
factor aS — (3'j on the right-hand side, i.e. by the non-zero determinant of 
the transformation between {A, B) and {A' , B') and so it is in fact the same 
equation. This shows assertion 1 of the theorem. 

Assume now that the equation ( |4.5D is associated with a pair {A,B). Con- 
sider the vector X = {Xi,X2,X^y G Y{? of the coefficients of at the 
highest derivatives Krr, Krw, K^^u respectively. Then 

Xi = A12B22 — A22B12 

X2 = A22B11 - A11B22 (4.12) 
= M1B12 - A12B11 

or 

X = AxB (4.13) 

where A = (An, yli2, A22)* and B = [Bu, B12, B22Y are three-dimensional 
vectors depending on r and w. Hence, for a fixed (r, w) both vectors A and 
B are orthogonal to X. The coefficients at Kr, Kyj and K yield equations 
which are differential consequences of (|4.12| ) and so they do not impose any 
additional restrictions on {A,B). Suppose now that there exist matrices A' 
and B' satisfying the cyclic condition (|2.3| ) and associated with the same 
fundamental equation. This means that the equation ( [4.13| ) has another 
solution i.e. that X = A' x B' so that the vectors A' and B' are orthogonal 
to X and in consequence they are linear combinations of A and B: A' = 
aA -f- (5B, B' = 7^4 -|- 6B with some coefficients that may depend on r and 
w. For the corresponding matrices it immediately follows that 

A' = aA + I3B, B' = -fA + 6B. 

It remains to show that the coefficients a, (3, 7, 6 in fact do not depend on r 
nor w. This can be shown by inserting the explicit form (p^) and (^^) of 



entries of matrices A, B, A' and B' into ( 4.12 ). This shows assertion 2 of 



the theorem. QED 



V. Affine inequivalent forms of fundamental equa- 
tion 

In this section we are interested in characterizing all different types of two- 
dimensional qLN systems admitting two functionally independent integrals 
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of motion E and F which are quadratic in velocities, i.e. systems of the 
form 

q = M = -^A-^Vk = -^B-^Vl (5.1) 

where M is the force of the system. Every such system is described by a pair 
of matrices A{q), B{q) satisfying the cychc conditions (|2.3D and by a pair of 
functions k{q), l{q) satisfying A~^V/c = B~^Vl. We remind the reader that 
the functions k/ det{A) and 1/ det{B) satisfy the same fundamental equation 
with the coefficients completely determined by the matrix elements of A and 
B. 

Let us first consider how the qLN system = 6~^E = 2A{q + ^^l^-'^VA;) 
transforms under the affine transformation of coordinates 

q = SQ + h, S £ GL{2, R), /i G (5.2) 

where Q = {Qi, ■ ■ ■ , QnY- It is easy to see, that under the affine transfor- 
mation (O) the generating function E transforms as 



Eiq{Q), qm = Q'S'A{q{Q))SQ + /c(g(Q)), (5.3) 

where q{Q) = SQ + h and so q{Q) = SQ. It can be shown by a direct 
verification, that the transformed matrix 

Aq{Q) = S'A{qiQ))S (5.4) 



in (5^) also satisfies the cyclic conditions ( p7^ ) and therefore ( |5.3D generates 
a qLN system. This means that the qLN system = 2A{q + ^A~^Vk) is 
indeed invariant with respect to the affine change of coordinates ( |5.2| ) . 

Let us now consider the system (^.1|). Using ([5.4]) one can prove that the 
fundamental equation associated with the pair (A, B) of matrices is also 
invariant with respect to the affine transfromations (|5.2| ). This means, that 
we can simplify this fundamental equation by performing an appropriate 



affine change of coordinates. But Theorem |4.7| makes it possible to classify 
fundamental equations, and therefore the corresponding qLN systems, by 
classifying pairs of matrices {A, B). Instead of working with the coefficients 
of the fundamental equation we can thus work with linear spans {\A + /ii?} 
of A and B. Since the affine transformations do not change the polynomial 
degree of matrices A^ B, the set of all linear spans of A and B can be divided 
into affine inequivalent classes corresponding to different polynomial degree 
of A and B. Each equivalence class will be represented by the algebraically 
simplest pair of matrices obtained by the use of affine transformations and 
linear combining of matrices (since the latter leave the fundamental equation 
unchanged, see above). 

In order to be more precise we shall introduce some notation. By A^'^^ 
(i = 0, 1, 2) we will denote all matrices A which satisfy the cyclic conditions 
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(2.3) and have the highest degree of polynomial entries equal to i. So, for 
example, the general form of matrices in the class A^^^ is 



bqi + a —\hq 



- 

cqi + 7 



+ i/3 



with arbitrary constants (parameters) b, c, a, (3, 7. We will use the symbol * 
to denote matrix elements determined by the symmetry of a given matrix. 
Moreover, by [A^^\B^^^ = 0,1,2) we will denote the class of (non- 
ordered) pairs {A, B) of linearly independent matrices A, B such that one 
of the matrices belongs to j4^*^ and the other to i?*^-'-*. We have, of course, 
[A^*), B^^^\ = [A^^\B^'^'^] and so we have precisely six such classes. Obviously, 
all classes [A^'^\ B^^^\ are invariant with respect to the affine transformations 
( |5^ . Notice that if {A,B) £ [A(2),5(2)] we can kill the coefficient 02 at 
the second degree monomials in B by subtraction B B — ^A so every 
element of this class can be reduced to an element in [A^'^\ B^^^]. Thus 
we have to consider only five classes. It is easy to realize that the five 
classes IA^^\B^^^] with i > j, j < 2 are invariant with respect to the affine 
transformations ( |5.2[ ) and with respect to taking linear combinations of pairs 
{A, B). These two operations can now be used to find for every class a simple 
representing pair [A, B) which has a minimal number of free parameters (a 
simple representative). Consider for example the class [A^'^^ S^'^^]. The 
general form of matrices belonging to this class is 



A 



aiql + 61^2 + ai 
* 



B 



aiqiq2 - ^biqi - \c1q2 + (3i/2 
aiql + ciqi + 71 



(5.5) 



"2 hi 2 

72 



Translation by the vector h = ~'^{ci-,biY kills bi and ci in the matrix 
A. Since translations obviously preserve the form of B the above pair of 
matrices attains the form 



A 



aiql+ai -aiqiq2 + (3i/2 
* aiql + 71 



B 



a2 132/2 
P2/2 72 



with some new constants denoted by the same letters as in (|5.5| ). Further, 
the transformation A\-^ A — ^B kills the coefficient /?i in A. In case when 
/32 = we can still kill f3i in A by the linear transformation q = SQ with 



5 



1 



-t7i-/3i/2 

1 t 



where t G R must be chosen so that ai + tPi/2 / and det(S') 7^ which 
can be always done. Finally, we can divide both matrices by oi and 2q2 
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respectively. So a simple representative of the class [A^'^\ B^'^^] has the form 



A 



9i +71 



B 



1/2 (52/2 
^2/2 72 



(5.6) 



with four essential parameters. Both separable and driven systems belong 
to this class since B = for separable systems and B = diag(l/2,0) for 
driven systems. 

We perform a similar reduction for each class [A^^\ B^^^], i > j, j < 2. The 
results are presented below. 

It is also easy to see that one can pass from one invariant class [A^^^B^^'^] 
to another by specifying values of free parameters. For example, by setting 
ai = we obtain [A^^\ B^^'^] from [A^'^\ B^^^]; by setting 62 = C2 = we 
get i?^'''*] and so on as shown in Figure |l[ This figure presents — for 

all classes [j4(*\s(-^)] — complete results of simplification of a generic pair 
{A, B) belonging to each class with the use of linear combinations and affine 
transformations . 



Below we list the form of the fundamental equation (|4.q ) corresponding to 
the simple representative pair (A, B) of each class as given in Figure |^. We 
use the notation Ki = dK/dqi, Kij = d'^K/dqjdqi. 

a) For [.4(2), 



9192) 



= Kii (-71/32 + 71^2 - 27291^2 - /32g'l - 9192) 

+2K12 (a27i - «172 - aiqi - 72^2 + "291 
+K22 (ai/32 - ai92 + 2a29i92 + hql - ll) 

+3Ki (-/3291 - 27292 - 29192) + 3^2 (-ai + 2a29i + /3292 - 29!) 
-Qq2K 

b) For [^(2)^^(0)] 

= i^ii (-71/32 - 2729192 - /329i) + 2i^i2 (71 - "172 - 7292 + 9i) 
+K22 (ai/32 + 2qiq2 + l32q^) + 3Ki (-/329i - 27292) 
+37^2 (291 + (32w) 



(5.7) 



2Kii (72/3i - 71/32 + (-72 + Pi)qi + 92 - 9i) 



c) For 

= 

+4i^i2 (0271 - ai72 - aiqi - 7292 - 9192) 
+2^:22 (ai/32 - "2/31 + 0291 + (/?2 - «i)92 - 
+3i^i (-272 + /3i - 3gi) + 3i^2 (-2ai + /32 



9i) 



392 



6K 
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^(2) R(l) 



+ ai -qiq2 
-qm Qi + 71 



"2 -^2 + 02/2 

* gi + 72 



ai = 



62 = C2 = 



"^(2)^5(0)- 



g| + ai -9192 
-Qm ql + 71 



1/2 152/2 
fi2/2 72 



ai = 



q2 + ai -\qi + \l3i 
* 71 



02 -\q2 + 5/^2 

* 91+72 



&2 = C2 = 



92 -^91 
-igi 



1/2 P2/2 
(^2/2 72 



61 = ci = 







[1 1 [1/2 (52/2 1 
[0 71 J ' [ /32/2 72 J 




A(o),i?(o)' 







Figure 1: Classification diagram 
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d) For 



= 27291 i^ii + 47292^^12 - 2(gi/2 + /?2g2)i^22 + 672^^1 - (5.8) 
e) For 

= 7i/32-fS^ii + 2(ai72 - a27i)-?^i2 - ai/52-^22 (5.9) 

What we present here is an illustrative characterization of different types 
of fundamental equations in terms of matrix pairs {A,B). This provides a 
good intuitive description of the world of qLN equations and helps to specify 
where two particular classes — separable potentials and driven systems — 
belong. An alternative way of classifying qLN equations with two quadratic 
integrals of motion is to simplify the fundamental equation ( |4.5| ) with the 
use of affine transformations as has been done for the Bertrand-Darboux 
equation |^] (see Example |5.1| below). This may amount to a similar picture 
as we have presented above but the principles of simplification of the third- 
order polynomials at Krr,Krw and K^jyij are more difficult to discern. This 
is yet to be done. 

Example 5.1 The classification of types of the Bertrand-Darboux equa- 
tion with respect to Euclidean transformations leads to three forms of this 
equation which are separable in either elliptic, parabolic or Cartesian coor- 
dinates. According to Corollary |3.1| , if a potential two-dimensional Newton 
system 

dV .. dV 

oqi oq2 

with V = V{qi,q2) possesses a second integral of motion of the form E = 
q^Aq + k{q) with A G A^'^\ then it has the qLN form q = -^A~^Vk = 
— ^B^^W with B = where I is 2 x 2 identity matrix. Moreover the 
potential V must satisfy the Bertrand-Darboux equation (^). This system 
belongs to the class [^(2)^^(o)j ^j^-j^ /?2 = and with a2 = 72 = 1/2. The 
corresponding fundamental equation is exactly the Bertrand-Darboux equa- 
tion since in this case K = V/ det(i?) = 41/. The simplification procedure 
described above does not alter the form of the matrix B = and so the 
corresponding fundamental equation ( ^.T] ) attains the form 

= (^22 - 141)^1^2 + Vuiqf - + 71 - ai) - 692^1 + 691^2 (5.10) 

which has only one essential parameter 71 — ai. This form of the Bertrand- 
Darboux equation separates in the elliptic coordinates Q 
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(where ^m^ = ai — 71) which are its characteristic coordinates. In these 
coordinates the Bertrand-Darboux equation ( 5.1C| ) takes the form 



2ry 2^ 

and its general solution is 

f{0-9{ri) 



^2 _ ^: 



,2 



with arbitrary functions / and g. The specification oi = reduces the 
class [A^'^\ B'^'^^\ to [^(^) , The corresponding fundamental equation 

after the simplification procedure attains the form (5.8). In the case when 



B = the final form of B (after simplification) will be exactly the same 
(i.e. with P2 = 0, Q!2 = 72 = 1/2) and so the fundamental equation ( |5.8| ) 
reads 

= qi{Vu - V22) + 2q2Vi2 + 3Vi 

It does not contain any parameters now. This equation separates in the 
parabolic coordinates [§] = (,1], q2 = ^{^"^ — ff) in which it takes the form 

= ^. + ;^^ + ;^^. 

and has the solution 

with arbitrary / and g. Further specification 61 = ci = leads to the class 
[A^'^\ B^^'^\ of constant symmetric matrices. The corresponding fundamental 
equation in the course of simplification attains the form ( ^.S| ) which in the 
case B = (again, this form of B survives the simplification procedure - in 
this case just the diagonalization of A by a rotation) yields V^,; = 0. This is 
the case of the Bertrand-Darboux equation separable in (rotated) Cartesian 
coordinates. 

The well known Bertrand-Darboux theorem (known also as the Whittaker 
theorem) ^] says that potential Newton equations admitting a second 
integral of motion quadratic in momenta admit separation in one of the four 
coordinate systems: elliptic, parabolic, polar or Cartesian. The remaining 
polar coordinates do not belong to our scheme, since radially symmetric 
potentials V{qi + gf) in the function E = ^{qf + g^) + Viqf + qi) have the 
angular momentum J = q\q2 — (?2<?i as an integral of motion. This means 
that F = is the second (quadratic in velocities) integral of motion of 
every potential system with radially symmetric potential. But the function 
F = has l{q) = and therefore this case does not belong to our theory. 
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The above example indicates that the fundamental equation plays the same 
role in the theory of qLN equations as the Bertrand-Darboux equation does 
in the theory of separable potential forces M = —dV/dq. For separable 
potentials the characteristic coordinates of the Bertrand-Darboux equation 
determine the coordinates of separation which makes it possible to solve the 
corresponding Newton equations by quadratures. In Section 7 we prove a 
similar result for the class of two-dimensional driven systems by showing 
that the characteristic coordinates of the fundamental equation associated 
with a given driven system separate this system i.e. that in these coordinates 
it is possible to integrate the system by quadratures. The question whether 
the characteristic coordinates of the fundamental equation separate general 
qLN systems admitting two integrals of motion remains to be investigated. 
We have here to do with a much broader theory depending on five essential 
parameters while the Bertrand-Darboux equation depends on one parameter 
only. 



VI. Hamiltonian structures and complete integra- 
bility 



In this section we will establish a Hamiltonian formulation of two-dimensional 
qLN systems and discuss their complete integrability. Let us consider first 
the qLN system = = 2A{q + ^A~^'Vk{q)) generated by the function 
E = q'^A(q)q + k{q), q = (<?i,<?2)* with the 2x2 matrix A(q) satisfying the 
cyclic conditions ( [2.31 ). This system usually does not have any Lagrangian 
formulation and thus it does not have the standard Hamiltonian formula- 
tion. However we can always embed this system in a Hamiltonian qLN 
system in the five-dimensional phase space of variables {qi,q2,PiiP2,d) as 
the following theorem states. 



Theorem 6.1 (Hamiltonian form of qLN systems) Let 

1 



(6.1) 



= q+ -A'\q)V {k{q) + dXdet{A{q))) 

with q = {qi,q2Y be the qLN system generated by 

E = q*A{q)q + k{q) + dX det{A{q)) = E + dX det{A) 

with some constant X and with d £ H. Let also A4 be the extended 5- 
dimensional phase space of variables {qi,q2,Pi,P2,d) with pi = (ji, i = 1,2. 



Then the system {6A) is equivalent to 



Q 







-|G(g) 


P 


P 




|G*(g) 


-lF{q,p) 


M{q,d) 


d 




— p* 


-M\q,d) 






(6.2) 
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where Vm = {d/dqi,d/dq2,d/dpi,d/dp2,d/ddY is the gradient operator in 
M and where the 2x2 matrices G and F and the vector M are given by 



1 /dAoo dA 



A22 —A12 
-A12 All 



1 (dA22 dAii \ , 

Miq,d) = M{q) - ^dXA'^V {det{A)) 
with M{q) = —^A~^'Vk being the force of the qLN system = 6^E. More- 



over, the antisymmetric matrix Ha is Poisson and so (6.i) is the Hamilto 



nian formulation of (6.1 ). 



Notice that the matrix G obtained above is symmetric due to the symmetry 
of A. 



Proof: Since Vvid = (0, 0, 0, 0, 1)* the equation ( p^ ) yields q = p, p = M = 
— ^A'^V {k + dXdet{A)), d = 0, i.e. it reproduces (|6.1D. The matrix 11^ is 
antisymmetric and it is straightforward to verify that it satisfies the Jacobi 
identity in the phase space M. QED 

We remind the reader that the operator 11 : T*M — > TM mapping fiberwise 
the cotangent bundle T*M of M into the tangent budle TM is Poisson if 
the bilinear mapping {•, -jn : G°°{M) x G°°{M) G°°{M) defined for any 
pair of functions /, g : ^ R by 

{/,9}n= <VA4/,nVAi5> 

(where < •,• > is the dual map between cotangent and tangent spaces of 
M.) is a. Poisson bracket. 

Remark 6.2 In the hyperplane d = the solutions of ( |6.2| ) coincide with 
the solutions of q = — ^j4~^V/c. Thus our original qLN system q = — ^^~^VA; 



is in a natural way embedded in the Hamiltonian system (6.2). 



Proposition 6.3 The function 

E = p^A{q)p + k{q) + dX det(^(g)) 

is a Casimir function for the Poisson operator in l\6.S^ , that is Ha^m^ 
0. 

One can check this proposition by a direct verification. 



A statement converse to the second statement of Theorem 6.1 also holds. 
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Theorem 6.4 Let the antisymmetric matrix 








-lG{q) 


V 


n = 




-^F{q,p) 


M{q,d) 




— 


-M\q,d) 






(6.3) 



be a Poisson operator in the space of variables {q,p,d). Then 
1. G{q) must have the form 

G{q): 

(thus it is symmetric) with some constants a,b,c,a, (3,'j and so 



aqj + cgi + 7 aqiq2 + |gi + §92 - f 
* aq2 +bq2 + a 



G 



A22 
-A21 



-A12 
All 



(6.4) 



(6.5) 



for some symmetric matrix A{q) satisfying the cyclic conditions (2.5). 
In other words, H = Ha with Ha defined in l \6.^ ) and with 



A 



2. F{q,p) must have the form 



G22 —G12 
-G21 Gil 



3. M{q, d) must have the form 

M{q,d) = M{q) + dXN{q) 



(6.6) 



where —2AM{q) = Vk for some function k{q), so i/ det(G) 7^ then 
M{q) = -\A-^Vk, and where N{q) = -iA-iV(det(A)). 



Proof: The conditions {{qi,qj}n,qk}n + cycl = and {{qi,qj}ii,Pk]ii + 
cycl = (where "cycl" means the cychc permutation of expressions) hold 
identically due to the block structure of 11. The condition qj}n-, d}ji + 
cycl = yields the symmetry of G: G = G^ . Further 

A/ dG, dG.,\ X dM, 



= {{qi,Pj}u,d}u + cycl = -- [pi 



Let us denote the right hand side of the above equality by —^Rij. Notice 
that dMj/dd can not depend on d and so we have dMj/dd = XNj{q) for 
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some vector N{q) = {Ni{q), N2{q)Y which yields M{q,d) = M{q) + dXN{q) 
for some vector M{q). By taking hnear combinations of the conditions 
Rij = and using the symmetry of G and the antisymmetry of F we get 
the fohowing sets of equations 

dGii dG22 ^ ^G^^ dG: 



dq2 



dqi 

dGi2 _ 
dqi 

Fl2 



dGu 
dqi 



2Ni, 
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dq2 



2N2 



dG 



12 



dq2 

P2N1 -P1N2 



(6.7) 

(6.8) 
(6.9) 



The equations (6.7) show that Gu and A'^i depend only on qi and that 
G22 and N2 depend only on 52- The equations (|6.8| ) give dNi/dqi = 
d'^Gi2/dqidq2 = dN2/dq2 and so all terms in this expression must be equal 
to a constant a. Integration yields 



Ni = aqi + c/2, N2 = aq2 + 6/2 



5.10) 



where h and c are integration constants. Substituting ( |6.10|) into ( |6.7| ) and 
(6.8) and integrating we get (|6.4| ). If we now introduce the symmetric matrix 
A by the equality (|6.5|) and use (|6.7]) then (|6.9|) will attain the form (|6^). 



It is straightforward to check that with the above forms of F and G the 
conditions {{pi,Pj}n,Pk}u + cycl = and {{qi,Pj}ii,Pk]u + cycl = are 
satisfied identically. 

Further, the condition {{pi,P2}ii^ d}n + cycl = after some calculations 
attains the form 







dqi 



(G11M2 - G21M1) 



d_ 

dq2 



(G22M1 - G12M2) 



which means that in the vector 



G22M1 - G12M2 




^22 


— G12 




' Ml ' 


-G21M1 + G11M2 




— G21 


Gil 




M2 



AM 



-\Vk or M 



the mixed derivatives of its components are equal and so this vector is equal 
to the gradient of some function —\k{q)^ that is AM 
-\A~^Vk. 

Finally, by direct calculation we verify that N 
statement 3 of the theorem is proved. 



-\A~'^V{d&t{A)) and so 
QED 



Remark 6.5 This theorem generalizes the result of In particular, if we 
assume M = — Vy(g) then we recover the known second Poisson operator 
for separable potential systems |^]. 
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Notice that M is the force of the 2-dimensional qLN system (|6.1|). This 
means that every Poisson operator of the form (|6.3| ) is the Poisson operator 
for some qLN system of the form ( |6.]j ) . 

We are now in position to investigate complete integrabihty of qLN systems 
admitting two quadratic, functionahy independent integrals of motion. No- 



tice first, that Theorem 6.4 provides us with an alternative way of charac- 
terizing qLN systems generated by a quadratic integral of motion E: by 
starting with a Poisson operator of the form ( |6.3| ) we arrive at qLN systems 
generated by the Hamiltonian H{q,p,d) = d which admit a quadratic inte- 
gral E. In a similar way the following theorem characterizes all qLN systems 
admitting two independent quadratic integrals E, F. 



Theorem 6.6 (Poisson pencil) Consider the antisymmetric operator 










P 








M{q) + d\N^{q) 




— p* 


-M\q) - dXNj,{q) 






(6.11) 



where 



G, 



A22 
-A21 



-A12 
An 



B22 
-B21 



-B12 
Bii 



Ga - fJ-Gs, 



with both matrices A and B satisfying the cyclic conditions (2.5), 

fJ-Bii 



rj,! 1 fd{A22-fiB22) d{An 

[^Ii\l2 = - I — P2 



(with F 



dqi dq2 
F\ Fa = -FX, Fb = -Fy and 



-Pi 



[F. 



A\12 



f^[FB[ 



12 



— A'^V (det(A)) + -fiB-^V (det(B)) = Na - ^Nb- 



Then n„ is Poisson if and only if 



M{q) 



-A-^Vk 
2 



-B-^Vl 
2 



(6.12) 



for some functions k{q) and l{q)- Moreover, if we let 



n„ = Hi 





Xr<t 

2^A 


-iGa 


P 

M + dXNA 






Xr<t 

2^B 






dXNB 


— p* 


-M* - dXN^ 










-dXN% 






then both operators Hi and JI2 cifc Poisson and so 
Poisson pencil. 



Hi — IJII2 is a 
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Proof: According to the proof of Theorem |6.4| the matrix 11^ satisfies all 
the Jacobi identities except possibly for {{pi,P2}n^) rf}n^ + cycl = 0, since 
differs from Ua-^b = — fJ-^B by the form of M{q) only. Like in the 



proof of Theorem |6^ we find that {{pi,P2}n^i t^ln^ + eye/ = yields that 
the mixed derivatives of the components of the vector —2{A — fj,B)M are 
equal and so —2{A — ^B)M = V{k — fj,l) for some functions k{q) and l{q). 
By comparing coefficients at different powers of fj, we get —2AM = Vk and 
-2BM = VI and thus M = -^A-'^Vk = -iB'^Vl. 



Further, Hi = 11^ in the notation of Theorem 6.4 so it is Poisson. Easy 
calculation shows that 112 is Poisson too. QED 

The above theorem states that if M{q) is the force of a qLN system admit- 
ting two functionally independent integrals of motion then the matrix 11^ is 
a Poisson pencil. We will establish its Casimir function, which will be a poly- 
nomial in ^. This will lead to a bi-Hamiltonian chain containing the qLN 
system ( |6.1D . We will prove that this chain is completely integrable. In this 
way we will show that our original qLN system q = — ^yl^^V/c = —^B^^Vl 
can be naturally embedded in a completely integrable bi-Hamiltonian sys- 
tem. 



Proposition 6.7 Suppose that is Poisson, i.e. that ( G.li. ) is satisfied. 
Then the function 



Hf, = p\A - ij,B)p + k- fil + dX det{A - fiB) 
is a Casimir function for n„, i.e. TluVHu = 0. 



(6.13) 



Proof: This proposition is a consequence of Proposition |6.3| . If we modify 
the matrix 11^ by substituting the matrix ^ by ^ — iJ,B and substituting k 
with k — fil we obtain the matrix 



(6.14) 






2-^M 


p 

Mm 


— 


-M* 






where 



1 



{A - fxB)-'V{k - fil) - -dX{A - fiB)-'V {det{A - fiB)) . 



Due to Proposition 6^ the function (|6.13| ) is the Casimir of (|6.14D . But 
( |6.14| ) is in fact equal to 11^ since it can be verified that —^{A — fiB)~^\7{k — 
fil) = -\A-^Vk = -\B~^Vl and that {A - fxBy^V {det{A - iiB)) = 
A-^V (det(^)) - fiB-^V (det(5)). QED 
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Let us collect terms in at different powers of /u 

= p^Ap + k + dX det{A) + ^ (-p*Bp - / - dXY^ + /i^ {dX det{B)) 
= E + HF + fi^H 

with Y = B11A22 + B22A11 — 2i?i2j4i2. Then the above proposition gives: 
= U^VH^ = (Hi - fiU2)V (^E + fiF + i?H) = 

By equating to zero the coefficients at different powers of ^ we obtain the 
following bi-Hamiltonian chain: 



liiVE = 
HiVF = 112'^E 

= HaVF 



(6.15) 



Theorem 6.8 The bi-Hamiltonian chain ( 6.13^ is completely integrable, 



i.e. both non-trivial bi-Hamiltonian vector fields 

Vi = HiVF = HzV-E, V2 = liiVH = HzVF 



in ( 6.15 ) are completely integrable. 



Proof: (modification of the proof of Liouville- Arnold theorem [^) Con- 
sider the 2-dimensional manifold M = {x £ M : E{x) = Eq,F{x) = 
Fq,H{x) = Hq} in M. Poisson brackets of all pairs of E,F,H induced 
by both structures Ho and Hi are equal to zero, since the functions E, F, H 
all belong to the same bi-Hamiltonian chain. For instance {F, H}ui = < 
VF, HiV^ > = < VF, U2VF > = {F, F}u2 = with the second equal- 
ity being a consequence of the bi-Hamiltonian structure of V2. It follows 
that the Lie bracket [Vi, V2] of both vector fields Vi and V2 is equal to zero, 
[Vi,V2] = [HiVF.HiV^] = 0, since the mappings HjV {i = 1,2) are Lie 
algebra homomorphisms between the Lie algebra of vector fields on and 
the Lie algebra of all smooth functions on M with the Lie bracket defined 
by [/i,/2] = {/i,/2}n,. Moreover < VE,Vi > = < VE,U2VE > = 
{E, E}u2 = and similarly < VF, Vi > = < VH, V2 > =0 which 
proves that Vi is tangent to J\f. In the same way one can show that V2 is 
also tangent to M. Direct verification shows that Vi and V2 are linearly 
independent. We thus have a 2-dimensional sub-manifold M in M equipped 
with a pair of linearly independent, commuting vector fields Vi and V2. We 
can now apply the construction of Liouville- Arnold |^ and conclude that 
both Vi and V2 are completely integrable. QED 
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Corollary 6.9 The qLN system q = M 



-B-^Vl with two 



linearly independent matrices A and B satisfying the cyclic conditions { 2^ ) 
is completely integrahle in the sense that the trajectories of the system 



(6.16) 



coincide on the hyperplane d = with the trajectories of the completely 
integrahle 5-dimensional system 



q 




p 


p 




M 



q 
p 

d 



V2 = UiVH = U2VF. 



(6.17) 



Proof: Consider the vector field V2 from ( |6.15| ). Obviously 

V2 = UiVH = X det{B)UiVd + XdUiV (det(fi)) 
and so in the hyperplane d = we have 



2\d=0 



Xdet{B) 



P 
M 





(6.18) 



which means that the hyperplane d = is invariant with respect to the action 
of the vector field V2. The formula ( |6.18D also shows that in the hyperplane 
d = the vector field of the system ( |6.17| ) is parallel to the vector field of 
the system (|6.16| ) and so their trajectories must coincide. QED 



Thus we have shown that the system ( |6.16 ) is embedded in the completely 
integrable bi-Hamiltonian system ( |6.17| ). The trajectories of ( 6.16| ) stay on 
the intersection of invariant manifolds for (f 
Also, since we can now solve the system 



7[T\ ) with the hyperplane d = 0. 
|6.17| ) by quadratures the time 
evolution of the coefficient Adet(-B) in ( 6.18| ) can be calculated which makes 
it possible to solve the system (|6.16| ) by quadratures too. 



VII. New types of separation variables for driven 
qLN systems 

In this section we study an important class of two-dimensional qLN equa- 
tions called driven systems. We find for all such systems their separation 
variables and prove their integrability by quadratures. The variables of 
separation are of a completely new type: they consist of families of conies 
which are non-confocal in contrast with the classical separability theory for 
potential systems. 
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We remind the reader that we call a two-dimensional Newton system driven 
if one of the two differential equations depends on one variable only. By 
renaming the variables if necessary, we can always arrange for such a system 
to take the form 

qi = Mi{qi,q2), , . 

q2 = M2{q2). ^ ' 

The second equation can be solved on its own and its solution q2{x) then 
determines the equation for qi , which explains the name "driven" . A driven 
system always has one integral of motion F = ^1/2 — / M2 dq2, obtained by 
integrating the second equation once, but in general there need not exist 
any others. 



Here we shall consider driven systems that admit a quasi-Lagrangian for- 
mulation q = —■^A^^\7k{q). Here, as usual, A{q) is a non-degenerate 2x2 
matrix satisfying the cyclic conditions ( |2.3D , i.e., a matrix of the form 



A 



aq2 + bq2 + a -aqiq2 - |gi - §^2 + f 

-aqiq2 - |gi - §92 + f ml + cgi + 7 



(7.2) 



Such a system always has two functionally independent integrals of motion 
E = q^Aq + k{q) and F = ql/2 - J M2 dq2. 

By examining the second component of the equation q = — ^A^^VA;(g), we 
immediately see that a qLN system is driven iff 

Ai2dik - Aiid2k = 2det(^)M2(g2), (7.3) 

for some function M2(q'2) depending on 52 only. We can produce driven 
qLN systems with any given M2{q2) and A{q) by solving for k{q) in this 
equation. The case An = is degenerate and will be treated separately 
later (see Remark |7.6| ), so we assume from now on that An 7^ 0. 

We start by introducing separation variables for (^) as characteristic co- 
ordinates for (17. c 



Definition 7.1 Define curvilinear coordinates {u,v) = {u{q),v{q)) as fol- 
lows. Let n be a parameter indexing the family of characteristic curves of 
O given by 

Auiqix)) 
-Au{q{x)) 

and let v = q2- 



(7.4) 



In other words, the curves given by ( |7.4| ) are the coordinate curves of con- 
stant u. For a given matrix A these curves can be explicitely calculated. 



In Theorem 7.7 we will describe these curves more explicitly. Let us just 
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note for the moment that they are not parallel to the curves of constant 
V, because of the assumption An / 0. Thus the above description really 
defines a coordinate system (at least locally). There is some freedom in the 
choice of u, but this will not affect our results. By abuse of notation we 
will write f{qi, 92) and f{u, v) for the same function / expressed in different 
coordinate systems. 



Lemma 7.2 The general solution of ( [7.5| j is 

k{u,v)= f{u) + D{u,v)g{v), (7.5) 
where f is an arbitrary function, D = det{A) and 

5(92) = . f ■, [ M2{q2)dq2. 
^llW2j J 

Proof: Along each characteristic curve q{x) given by (|7.4|) we can consider 



Q as an ODE 

^Hqix)) = 2D{q{x))M2{q2{x)), 

with general solution 

k{q{x))=D{q{x))g{q2{x))+f, 

where / is a constant of integration. This can be verified by direct differen- 
tiation; the cyclic conditions imply that 

^D{q{x)) = diD qi + 820 q2 = diD Au + 820 (-An) = -D Ss^n, 
and thus 

^%) = l{Diq)g(q2) + f)='^9 + D02gq2 

= -Dd2Ang + Dd2g{-An) 
= -Dd2{Aiig) = 2DM2. 

The constant of integration / can be different for different characteristic 
curves, so when we express the result in terms of u and v, f will depend on 
u (but not on v). QED 



Lemma 7.3 Equation (7.S) is equivalent, under the substitution k = K det{A), 
to the equation 

AuduK - Andi2K - ^ SsAn diK = 0, (7.6) 



which is the fundamental equation ( 1^.5[ ) associated with the matrices A and 


1/2 
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Proof: Equation { ^T^ implies 



^ f -Audik + And2k\ ^ 
[ d^) ) = 

Conversely, this expression can be integrated to give ( |7.3| ), where M2(g2) 
is an arbitrary function of integration. By substituting k = K det{A) and 
simplifying the resulting expression using the cyclic conditions one obtains 



(7.6). Comparison with the general expression for the fundamental equation 



in Theorem 4.1 proves the second statement of the lemma. QED 



Remark 7.4 The fundamental equation ( [7.6| ) is hyperbolic. Its character- 
istic coordinates are precisely the coordinates {u, v) of Definition The 
general solution is K{u,v) = f{u)/D{u,v) + g{v), as can be seen by com- 
bining the above lemmas. 



Let us turn to the question of how to integrate a driven qLN system. The 
solution q2{x) of the second equation can be found by quadrature from 
F = q^/2-jM2dq2: 

dq2 



dx. (7.7) 

/2F + 2 / M2 dq2 
Inserting q2{x) and q2{x) into 

E = Aii{q2)ql + 2Ai2{qi,q2) qih + ^22(^1) + ^(91,^2) 

would give a first order ODE for gi(x), but there is no obvious way to solve 
this equation since the variables qi and x do not separate. We will now show 
how to proceed instead. 



Theorem 7.5 Every driven qLN system can he integrated by quadratures 
using the characteristic coordinates (u, v) of the fundamental equation ( |7. 6| j 
as separation variables. 



Proof: We use the notation of Lemma [7.2| . Let the system be generated 
by = q^Aq + k{q) with k{u,v) = f{u) + D{u,v)g{v). Since v = q2, we 
can express E as E = ^ (u^ + An{v)g{v)) and calculate v{x) by quadrature, 
as above. Now note that since the curves of constant u by definition have 
tangent q = (^12,-^11)*, we must have Vn = /3(g) (^n, A12)* for some 
function p{q), whose exact form depends on the choice of u. This gives 
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ii = diuqi + d2uq2 = p{q) {Auqi + ^12^2), and thus 



P^Au 



P^An 



P^An 



p'An 
p'A,, 



Axxix + 2^129192 + 

E - A22ql - k{q) + -^g 



E 

E- 



A22 - ^) - /(^) - D9{v) 



D 
A~ 



11 



(2F - Angiv)) - fiu) - Dg{v] 



E-^F-f{u] 
All 



In order to complete the proof, we will show that p{u, v) = 4>{u) \Aii{v)\~^^'^ 
and D{u,v) / Aii{v) = ip{u) for some functions cp and ■0, since this implies 
that the variables u and x separate. Explicitly, we can then find u{x) from 
the quadrature 



du 



^{u)^E-2^P{u)F - f{u) 



dx 



Aii{v{x)y 



(7.8) 



after which the inverse coordinate transformation gives us qi{x). Notice, 
that for a given matrix A the characteristic coordinates (n, v) can be cal- 
culated explicitly so that the function p and thus (p and V can be easily 
calculated and used in the quadrature (^) above. The theorem covers 
however all the cases at once without any need of calculating p explicitly. 

To see that p{u,v) = (j){u) 17411(^)1^^^^, note that d^u = d2iu implies that 
p{q) satisfies the PDE 



= di{pAi2) - d2{pA 



u) = Ai2dip - Aiid2P - -52^11 p, 



which has the same characteristic curves ( [7.4[ ) as equation ( |7.3| ) . Along such 
a curve we determine p by integrating 

d 3 

— = -d2Aii{q2{x)) p{q{x)) 
which, taking into account q2{x) = —Aii{q2{x)), gives 

p{qix)) = cP \Aiiiq2ix)r'/\ 

The integration constant (f) can be different on different characteristic curves, 
so changing to (n, v) coordinates we obtain 



p{u,v) = (j){u) \Aii{v)\ 



-3/2 



34 



as desired. 



Finally, we calculate the total derivative of the function ^p{q) = D{q) /An (^2) 
along a characteristic curve: 

= ft - ^) A. - ft (a.. - 1^) An. 

Using the cyclic conditions, we find that this expression is identically zero. 
This implies that tp is constant along the coordinate curves of constant u, 
i.e., ^/J = ip{u). This completes the proof. QED 



Remark 7.6 The degenerate case An = can be treated as follows. Since 
a = b = a = 0, the expression ([7^) for A reduces to 





§92 + 1 



--2Q2 + f 

cqi + 7 



(7.9) 



and equation ( [7.3^ reduces to 

Audik = 2(^-Al^)M2{q2), 

with the general solution 

k{q) = -2 Ai2{q2) M2{q2) qi + k2{q2). 

We calculate q2{x) by quadrature as before. Inserting q2{x) and q2ix) into 
E = 2Ai2qiq2 + ^22^2 + k{qi,q2) yields in this case an equation of the form 
^1(2;) + ^{x)qi{x) = r]{x), from which we can find qi{x) by quadrature. 



Theorem 7.7 The separation coordinates for driven qLN systems, i.e., the 
characteristic coordinates {u, v) of the fundamental equation ( [7. 6| ^, are of 
one of the following types, determined by the coefficients in the matrix A: 

1. Fan-like hyperbolic, if a ^ and 6^/4 — aa = 0. 

2. Axial hyperbolic, if a ^ and b'^/4 — aa < 0. 

3. Two-point elliptic-hyperbolic, if a ^ and 6^/4 — aa > 0. 

4. One-point parabolic, if a = and 6 7^ 0. 

5. Parallel parabolic, if a = and 6=0. 
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Proof: We will compute explicitly the curves given by (7.4), which consti- 
tute the curves of constant u. (The curves of constant v are just horizontal 
lines, since v = q2-) Inserting the explicit expression ( |7.2| ) for the matrix A 
into (|7.4|), we obtain 



^1 
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aqiq2 + |gi + f 92 - f 
aq^ + bq2 + a 



(7.10) 



When solving these equations, we distinguish four different cases, depending 
on the values of the parameters in A. 



The case a / 0. 

By setting ri = aqi + c/2 and r2 = aq2 + b/2, which is just rescaling of the 
axes and translation of the origin, we transform ( [7.101 ) into 



ri 




-rir2 + Ci 


r2 




-rl + C2 



where 



Ci = 6c/4 + a/3/2 
C2 = 6^4 - aa 



(7.11) 



Subcase C2 = (type [![). 

Either ra = 0, or ra = (x + Di)-^ and n = Ci{x + Di)/2 + D2{x + Di)-^ , 
where Di and D2 are constants of integration. Eliminating x and writing u 
instead of D2, we obtain 

ri = ^+nr2, (7.12) 
2r2 

which represents a family of hyperbolas, each with asymptotes ra = and 
r2 = ri/u. The solution ra = found above corresponds to the limiting 
cases u — > ±00 (see Figure 2). 



Subcase C2 / (types |^ and 

The substitution si = ri — Cira/Ca, S2 



r2 yields 



and thus 



resulting in 







-S1S2 


S2 




-sl + C2 


dsi 


_ '^1 


S2 


dS2 


S2 


S2-C2 



sl 



(7.13) 



(7.14) 



If C2 < (type §), this represents in the s-plane a family of hyperbo- 
las centered around the si-axis, with asymptotes S2 = isi/u and vertices 
{±u^/^C^,0) (Figure 3a). 
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If C2 > (type |3|) , we obtain in the region | S2 1 > \/C^ a family of hyperbolas 
with asymptotes S2 = ±si/ii and vertices (0, ib\/C2^), and in the region 
|s2| < VC2 a family of ellipses with vertices (0, ±v^C2^) and (0, zttx-y/C^) 
(Figure 4a). The corresponding curves in the r-plane are obtained by a 
shear in the si-direction with factor C1/C2 (Figures 3b and 4b). They are 
still hyperbolas and ellipses, but not aligned parallel with the r-axes. 



The case a = 0. 

Subcase 6/0 (type 

Translating the origin by ri 
obtain 



which yields 



qi — {ac/b'^ + P/b) and r2 







r2 





-br2 



ri - Tr2 
b 



ur2- 



q2 + a/b, we 
(7.15) 

(7.16) 



With si = ri — cr2/b,S2 = ^2, we obtain in the s-plane a family of parabo- 
las 52 = Si/u (Figure 5a). The corresponding curves in the r-plane are 
parabolas obtained by a shear in the si-direction with factor c/b (Figure 
5b). 

Subcase 6 = (type 

Here we can assume that a 7^ 0, or else we get the degenerate case An = 0. 
A simple calculation shows that 



c 2 /3 



(7.17) 



which is a family of translated parabolas seen on Figure 6 (or straight lines, 
if c = 0). 

QED 




Figure 2. Fan-like hyperbolic. 
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Figure 3b. Axial-hyperbolic in the r-plane. 
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Figure 5a. One-point parabolic in the s-plane. 




Figure 5b. One-point parabolic in the r-plane 



40 



Figure 6. Parallel-parabolic. 



VIII. Examples and applications 

The notion of a qLN force M{q) = —^A^^{q)^k{q) naturally generalizes 
the concept of a potential force M{q) = —'Vk{q), which is a special case. 
The qlN forces admit an integral of motion quadratic in velocities, which in 
the potential case becomes the energy integral. The function k{q) may be 
called a "quasi-potential" of the force M{q). 

A given force is easy to test for the existence of a qLN formulation, pro- 
vided that one knows the general form of the matrix A{q) solving the cyclic 
condition ( |2.3D . In two dimensions A{q), given by 

. , aw"^ + hw + a —arw — |r — + f 

A{r,w) = h r 8 9 ! 

—arw — ~ 2'^ + 2 ^ + cr + 7 

depends on 6 arbitrary parameters, and a qLN formulation exists provided 
that the mixed derivatives of V/c(r, u;) = —2A[r,w)M{r,w) are equal for 
some non-zero values of the parameters a, 6, c, a, /5, 7. We thus have the 
following lemma: 
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Lemma 8.1 A given force M{r,w) = {Mi{r,w), M2{r,w)Y admits a qLN 
formulation M{r,w) = —^A~'^{r,w)'Vk{r,w) if and only if there is a non- 
trivial solution A, with det(A) ^ 0, of the equation 

d ( bed 
= ("'^^ + 6w + a)Mx + {-arw - -r - -li; + -AM^ 

aw \ 2 2 2 

d f , b C P XT, r / 2 \,r 

— [[-arw - -r - -w + -^)^'h + [ar + cr + 7)Af2 

Lemma 8.2 (Criterion of integrability, n=2) Equation ^. j| j has a two- 
parameter family of solutions for A{r,w) if and only if q = M{q) admits two 
functionally independent integrals of motion E and E quadratic in velocities. 

Proof: If such E = q^Aq + k and E = q^Bq + I exist, then \E + /xF = 
g*(AA + iJ,B)q + {\k + ^l) \s an integral of motion for all A and /x, and thus 



\A + is a two- parameter solution of 8.1. 

Conversely, if there is a two-parameter solution D{\, /i) of 8.1, then there are 
linearly independent integrals E and E with A = D(l, 0) and B = D{0, 1). 
QED 

These two lemmas make it simple to test a given two-dimensional force 
for the existence of a qLN formulation, and to show integrability if a two- 
parameter family of solutions for A exists. 

Example 8.3 (gH-H system). The generalized Henon-Heiles (gH-H) sys- 
tem [|^] is defined by the potential 

T/r \ 2 C2gi Co , 
y{Qi,Q2) = ciqiq2 3" + 2g|' '^i'^2/^- 

It is known to be integrable in three cases: the Korteweg-de Vries (KdV) 
case 6ci + C2 = 0, the Sawada-Kotera (S-K) case ci + C2 = and the Kaup- 
Kupershmidt (K-K) case 16ci + C2 = 0. In the KdV case, and also in the 
S-K case if cq = 0, the second integral of motion is quadratic in velocities, in 
the other cases it is quartic. The system appears naturally when integrating 
the equation = {jd^ + 2cQud + cqUx){\uxx — jC2u'^), which for the above 
cases corresponds to the stationary flow of the fifth order KdV, S-K and K-K 



soliton equations. This observation explained |11] the remarkable connection 
of the integrable cases of the gH-H system with soliton hierarchies. 

We shall apply our criterion for existence of a qLN formulation to two New- 
ton representations of the gH-H system; the original system in (7- variables 

dV 2 2 

qi = - = -ciq2 +C2qi, 



92 = -T^ = -2ciqiq2 + ^, (8.2) 
dq2 q^ 
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and another system in r-variables 

r2 + (ci + C2)rl, 



r2 = C3- 10cirir2 - lOci [ci + — ]r{, 



C2 



i.3) 



which is equivalent IQ to the g-system under the map ri = qi, r2 = —ci{q(+ 



4) 



-4ci 



The r-system does not have any 



natural Lagrangian or Hamiltonian formulation and its integrability has 
previously been studied only through its equivalence with the original gH-H 
system. We will show here a more direct approach based on the qLN theory. 

Beginning with the g-system, we insert the right-hand side M from (8.2) 
into (pni), identifying (^1,92) with (r, it;) as usual. Since the powers r^w^ are 
linearly independent, the coefficients at different powers must be individually 
zero. This gives, after some simplification, that o = 6 = 0, a = 7 arbitrary, 
(6ci + C2)c = 0, (ci +02)/? = and co(3 = 0. This means that we always have 
a solution A = tl (of course corresponding to the energy integral, since the 
system has a potential). Moreover, in two cases there exists a two-parameter 
solution; when 6ci + C2 = 0, 



A = t 



and when ci + C2 = cq = 0, 



1 





+ s 



-q2 

-q2 2gi 



" 1 


" 




" 


1 " 






+ s 


1 








1 



A = t 



So in this way we have recovered the KdV and S-K (cq = 0) cases, while 
the K-K and S-K (cq / 0) cases, with a quartic extra integral, fall outside 
of the qLN theory. 

Performing the same procedure for the r-system (|8.3|), we find that a = 0, 



4ci7-|-6 = 0, (31ci-|-C2)c = {2>Ci+ciC2)c = 0, 2a — Scsc = and (6ci + C2)/3 



0. Since we have excluded the trivial case ci 
so that the solution is 



0, it follows that c = a = 0, 



-2r2 ri 
n l/2ci 



A = t 

except for the KdV case 6ci + C2 = which admits a two-parameter solution 



-2r2 ri 
n l/2ci 





' 


1 ' 


+ s 


1 






A = t 



This agrees with the known fact ||12|] that for the r-system the second integral 
is quartic in velocities in the S-K and K-K cases and thus cannot be found 
by this method. 
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Suppose, however, that we had found the second integral in these cases by 
some other method. Then we would still have use for the qLN theory in 
proving the system's integrability, since the r-paramerization admits the 
nonstandard Hamiltonian formulation (6^) with A = 1 and C3 playing the 
role of the fifth variable d. (This coincides with the Hamiltonian formulation 



that was found in |12] by transferring the standard Hamiltonian formulation 
from the g-parametrization, except for naming the momenta in reverse order; 
here pi = ri, while in that paper si = r2, -52 = ri-) 

For example, in the S-K case (ci = — C2 = 1/2) we have 



n 




" 


-1/2 


-ri/2 


Pi 









-ri/2 


r2 


P2 


Pi 




* 





-Pi/2 


r2 


P2 




* 


* 





C3 - 5rir2 - |rf 


C3 . 




* * 


* 


* 






yA/(C3, 



with the commuting integrals of motion E = — 2r2pf + 2ripip2+P2 + 4''i''2 + 
|rf(r^ _l_ 5^,^) + C3(— — 2r2), which is a Casimir, the Hamiltonian C3, and 
F = ^pf- Grir2Pi + (3rf - r2)piP2 + np^ + Irfr^ + W-rjr2 + |rf + ir| + 



(— 2rir2 — rf + ^Pi)cs, which is quartic in momenta. 



Example 8.4 (Quadratic force) As an example of direct application of 
the criterions in Lemma 8.1 and Lemma 8.2, let us determine when the 



simple homogeneous quadratic force system 



f = 6ir'^ + up' , 
u) = + 92vp 



admits two quadratic integrals of motion, so that it is integrable in the 
sense explained previously. Plugging into equation ( ^.Ij ) and equating the 
coefficients of independent powers r'^w^ to zero, we find a = 0, 59ib — 7c = 0, 
7b - 5020 = 0, -6I26 + Oic = 0, eip + 27 = 0, 2a + 6*2^ = 0. The two last 
equations always give rise to a one-parameter solution for a, /?, 7, while the 
equations involving b and c admit a (real-valued) one-parameter solution if 
and only if = 62 = 7/5, which thus is the only case admitting a two- 
parameter solution for the system as a whole. The solution is in this case 



A = t 



w 

-(r + w)/2 



-{r + w)/2 
r 





' -7 


5 ' 


+ s 




-7 


5 



Integrating Vfc = —2AM we find the two independent integrals of motions 
E = wr'^ — {r + w)rw + riiP + (5r'^ — 12r^w + lAr'^w'^ — 12rw'^ + 5w^)/28 
and F = -7r'^ + Wrw + -7w'^ + A{r^ + w^)/5. 
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The next example illustrates how an application of our recursion formula 
( 4.10| ) reproduces a known [13| family of potentials separable in parabolic 
coordinates. 



Example 8.5 (Parabolic separable potentials) The family of integrable 
potentials 



Vnir,w) 



with second integral /„ 




-wr 



f.2 _|_ j.'iy ^a^s found in 



. These 

potentials are all separable in parabolic coordinates r = ^rj, w = (^^ — 
r/^)/2 by the Bertrand-Darboux theorem and include V3 = Ar'^w + 8w^, 
which is the KdV case of the Henon-Heiles potential, as well as V4 = r'^ + 
12r'^w'^ + 16w^, the so-called "l:12:16-potential." We identify the matrices 



A 



—w r/2 
r/2 



B 



1/2 
1/2 



in the integrals of motion, and take ki = r Vq = r , li = Vi = 2w. It is then 
easy to show by induction that our recursion formula ( 4.1C| ), which in this 
case reads km = —r'^lm-i, Im = —2wlm-i — km-i, reproduces this family up 
to a sign: (-1)"'-'^ km = r^Vm-i, ' 



= r^Vm-i, i-ir-^lm = for all m > 1. 



In order to give an impression of the wealth of different types of non-potential 
Newton forces belonging to our theory, we will now examine solutions of the 
fundamental equation (|4.5| ) for some specified pairs of matrices A and B. 
Any such solution corresponds to an integrable qLN system, and once one 
solution has been found, a whole family of solutions can be constructed using 



the recursion theorem 4.3 



Example 8.6 (One-dimensional complex motion) If we take A and B 
as 



A 







-1 



B 



1 

1 



then the fundamental equation reduces to the Laplace equation Krr+Kw^i, = 
0. Given a solution K(r, w), i.e. a harmonic function, we have k = K dei{A) 
—K, so k is also harmonic. We find the corresponding / from the relation 



V/ = BA-^Vk 




1 



-1 







which is nothing but the Cauchy-Riemann equations for k and Z, so Z is 
the harmonic conjugate of k. The corresponding qLN system r = —kr/2, 
w = kw/2 can be integrated by introducing the complex variable z = r + iw 
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and the complex integral of motion £ = E + iF = [r"^ — uP' + k) + i{2rw + l) = 
{f + iw)'^ + {k + il) = + /(z), where /(z) = k{z) + il{z) is analytic. We can 
now determine z, and thus r and w, from i = ^\JE — f{z) by one complex 
quadrature. 

Repeated application of the recursion formula ( [4.10 ) yields in this case the 
standard cycle of conjugate harmonic pairs {k,l) {l,—k) {—k,—l) — > 
{-l,k)^{k,l). 



Example 8.7 (Fundamental equation separable in polar coordi- 
nates) Let 



—2w r 


, B = 


' 


w 








r 


w 


-2r 



Then the fundamental equation becomes = 2{r'^Krr + 2r'wKrw+'w'^Kwii,) + 
9{rKr + wKyj) + QK, which in polar coordinates (r = Rcoscp, w = Rsincp) 
transforms into = 2B?Krr + 9RKji + QK. The general solution of this 
equation is K{R, <j)) = M(t>)R''^ + go{(t>)R~^^'^ , for some arbitrary functions 
/o and go- Changing back to r and w, we find that the general solution of 
the fundamental equation in this case is 

K{r, w) = /o arctan — + go arctan — — , 

V rjl + {w/ry V '^Z (1 + (u>/r)2)3/2 

where / and g are arbitrary functions. 



Example 8.8 Let 

h 
62 

for some (non-zero) constants oi, 02, 61, 62- Then the fundamental equation 
becomes = —h2rwKrr + {hir"^ — b2W^ + 6201 — bia2)Krw + cirwKww — 
3b2wKr + SbirKuj. Finding the general solution seems difficult in this case, 
but looking for a particular solution of the form K{r, w) = f{r) + g{w), we 
find 

K{r, w) = co(6ir^ + 62^^^) + ^ H o + C3, 

where the Cj are arbitrary constants. To simplify formulas we let cq = 
1, ci = C2 = C3 = 0, which gives Ko{r,w) = bir"^ + b2w'^ . We take 
li = Kodei{B) = 6162-^ and determine ki from Vfci = AB^^Vli., which 



A 



-w^ + ai 
rw 



rw 

.2 1 



^2 



B 
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gives ki{r,w) = 6i&2(oi^^ + a2w'^). The corresponding qLN system q = 
-^A'^Vki = -\B~^Vli reads 

r = -bib2r, 
iv = —bib2W. 

This system is separated in r and w, and therefore trivially integrable. How- 
ever, we can obtain a more interesting integrable system by using the recur- 
sion formula (^^), which here gives 

0102 - oir^ - a2w'^ , 
6162 

= {aia2 - air"^ - a2w'^){bir'^ + b2w'^), 

, 0162 + a2^i - bir^ - b2uP' , , 

h = TT h - h 

&162 

= a2bjr'^ + aiblw^ - {bir^ + b2W^)'^ . 

The corresponding integrable qLN system q = —^A~^Vk2 = — ^-B^^V/2 is 

r = 26ir'^ -|- 2b2rw'^ — 026ir, 
id = 262^^ + 261 wr^ — aib2W. 

We derive from this also another solution Ki = I2/ det{B) = /2/&1&2 of the 
fundamental equation, not of the form Ki = f{r) + g{w). Having found 
Kq and Ki, we can employ the two-step recursion ( |4.11[ ), which in this case 
reads 

_ 0162 -I- 0261 - bir"^ - b2w'^ 0102 - oir^ - a2tf^ 

61 &2 O1O2 

to find a whole sequence {Km}°^^ of solutions. The first few are 
bib2{air^ + a2w'^) 



01O2 — oir^ — 02^2 ' 



= bir^ + b2W^, 



a2b\r'^ + oiftaU)^ - (6ir^ -|- 62^;^)^ 



6162 



-?^2 = {{bir^ + b2W^y + {a2b{r^ + ai¥2w ){aib2 + a2bi - bir^ - b2W 
-{a2blr'^ + aiblw"^ + aia2bib2){bir'^ + b2w'^)) /blbl. 



Finally, let us conclude with an example of a three-dimensional qLN system. 
A detailed treatment of higher-dimensional qLN systems is presented in a 
separate article Q. 
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Example 8.9 The Newton system 



ri = — lOr^ + 4r2, 

r2 = — 16rir2 + lOrf + 4r3, 

rs = — 20rir3 — 8r| + 30rxr2 — 15rf + d 



was found in as a parametrization of the 7th order stationary KdV flow. 
It has three integrals of motion 



y.2 

^ 2 2 2 S ^ 

£Ji = rir^ ~^ ~2 ~^ lOr^^rs — 4r2r3 + 8rir2 — 10r]^r2 + 3ri — dri, 
E2 = r^rf - rifl + r2rir2 - ?^2?^3 - ririr-^ + 

4r;^r2 + 5r]^r2 — — r^^ — 4r2 + — 12rir2r3 H — ^ — h ar2, 



1 



2 -2 , 2 -2 , -2 

rari + r^ra + r3 



(2rir2 + 4r3)rir2 + 2rif2r^ + 2r2fir-i) + 



rir| - 3rjr| + ^rfr2 + 2rir| + ^rfra - rlr2r^ + r^r^ - ^(rir2 - rs), 

which all are quadratic in velocities. This means that the system is generated 
by any of them through the quasi-Lagrangian equations. From the velocity- 
dependent parts we find 



" 





1 " 




" 2r3 


r2 


-n 




r| 


-rir2 - 2r3 


r2 





1 







r2 


-2ri 


-1 




-rir2 - 2r3 




n 


1 










-ri 


-1 







r2 


ri 


1 



as examples of 3 x 3-matrices satisfying the cyclic conditions (2.3). 



IX. Conclusions 

In this article we have developed a new theory — the theory of quasi- 
Lagrangian Newton equations. It was originally inspired by interesting prop- 
erties of the second stationary flow of the Harry Dym hierarchy, which led 
us to a broad theory which encompasses the classical separability theory 
but goes far beyond the classical results — the classical Bertrand-Darboux 
theory of separability for two-dimensional potential forces depends on one 
essential free parameter while our theory depends on five parameters. 

The main part of this work has been focused on two-dimensional qLN sys- 
tems which admit two integrals of motion E and F quadratic in veloci- 
ties. These systems have only a non-standard Hamiltonian formulation and 
are completely integrable by embedding into five-dimensional Liouville in- 
tegrable systems. All such qLN systems are characterized by a single PDE 
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called here the fundamental equation. We have shown that there is a one- 
to-one correspondence between fundamental equations and linear pencils 
-|- fj,B of matrices A and B. These linear pencils have been classified 
in Sec. 5. In Sec. 7 the class of driven systems has been studied in de- 
tail and new types of separation variables (non-confocal conies) have been 
found. We have also shown that any given force can be effectively tested for 
the existence of qLN formulation, which can further be used for unveiling 
its complete integrability and for solving the corresponding Newton equa- 
tion. We have illustrated this by several examples including the generalized 
Henon-Heiles system (Sec. 8). 

There are several natural directions of development of the theory of qLN sys- 
tems, n-dimensional versions of our main theorems on fundamental equation 
and on complete integrability have already been formulated and proved in 

§■ 

The great wealth of different types of integrable Newton equations contained 
in the fundamental equation remains to be studied. Here we have only 
discussed two special cases: separable systems and driven systems. However, 
one of the most challenging questions yet to be answered is the existence of 
separation variables for the fundamental equation in its most general form. 
It can lead to new and interesting connections with the classical theory of 
separability of the Hamilton-Jacobi equation and of linear PDE's. 
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